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Abstract. We present the results of a broad band (3-200 keV) spectral analysis of BeppoSA'K and XMM-Newton 
observations of the Vela plerion. The hard X-ray (>15 keV) emission is found to be extended over a region of 
12' — 15' radius, corresponding to a size of about 1.0-1.3 pc for a source distance of 290 pc. A single power law 
does not give an acceptable fit while a broken power law or a log-parabolic law nicely fit the data. The former 
spectral model has the photon index Pi = 1.66±0.01 for energies lower than the break value equal to 12.5 keV and 
the index F 2 = 2.01±0.05 up to about 200 keV. The total X-ray luminosity of the Vela plerion is Lx = 5.5 x 10®^ 
erg s“^, which implies a conversion factor of the spin-down power of ~10“^. 
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1. Introduction 

The structure and evolution of filled-centre SNRs, also 
known as plerions according the original definition by 
Weiler & Panagia (1978), is a rich subject that involves 
many physical processes from hydrodynamics to high en¬ 
ergy particle acceleration and radiation. Plerions are pow¬ 
ered by young pulsars which inject in the remnant a rela¬ 
tivistic wind containing high energy particles and a mag¬ 
netic field. The physical mechanisms that convert the orig¬ 
inal Pointing-dominated energy flux into particle energy 
is not completely understood. Furthermore, there is evi¬ 
dence that particles are accelerated inside the plerion at 
a relatively large distance from the pulsar but the basic 
mechanism is not entirely clear. Plerion emission is ob¬ 
served across the entire electromagnetic spectrum, from 
radio waves to the most energetic 7-rays. However, in some 
frequency ranges, only little data are available: for exam¬ 
ple, in the hard X/softy-rays very little is known about 
the spectra of plerions. 

The Vela SNR is one of the nearest plerions and its 
structure can be studied in detail. Its distance is esti¬ 
mated to be about 300 pc (Caraveo et al. 2001, Dodson 
et al. 2003), and therefore an angular size of 1' corre¬ 
sponds to 0.09 pc (2.7x10^^ cm). The best Chandra X- 
ray images show a compact nebula with a jet-counter jet 
originating from the pulsar and a double bow structure 
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(Helfand et al. 2001, Pavlov et al. 2001a). Kargaltsev & 
Pavlov (2003) presented maps of the 2-10 keV photon 
index {F{E^) oc E~^) with resolutions of 2".5 and 10" 
which show that the spectrum is harder (P ~ 1.2 — 1.4) 
in the near surroundings of the pulsar and becomes softer 
(F ~ 1.6—1.8) at distances of 2'—3'. A more detailed study 
of the inner Vela Pulsar Wind Nebula (PWN), based on 
13 X-ray images, spread over a time interval of 2.5 years, 
has been performed by Pavlov et al. (2003). In particular, 
these authors report the presence of bright blobs mov¬ 
ing away along the outer jet and fading on time scales of 
weeks. Photon indices in the band 1.0-8.0 keV derived for 
the various structures are generally between 1.2 and 1.5. 

Results of a spectral analysis of XMM data have been 
presented by Mori et al. (2004): their spectral analysis 
of PN data, extracted from a circular region around the 
pulsar position within a radius of 1', gave a photon index 
of 1.64±0.08. 

The present knowledge about the hard X-ray emission 
from the Vela SNR is rather poor. The most recent data 
are those obtained in the energy range 60-400 keV with 
OSSE-CGRO (FOV 3°.8xll°.4) (De Jager, Harding & 
Strickman 1996): a single power law best fit indicated a 
quite hard spectrum with the rather uncertain photon in¬ 
dex of 1.6±0.5. It was found consistent with the previous 
(2-25 keV) data obtained with the Birmingham Spacelab 
2 coded mask detector (Willmore et al. 1992), that gave 
F = 1.74±0.08 inside a circle of radius 6'—7'. Furthermore 
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Fig. 1. Left panel: XMM-MOSl image of the Vela PWN in the 3-10 keV band; contour levels are at 5x lx 10“^, 
2x 10“^ and 8x10“^ times the brightest point. Right panel: MEGS image of the source; outer contour levels are at 
3x 10“^, 7X10“^, 13x 20“^ and 33x 10“^ times the brightest point. Note the more extended emission in the south-west 
and south-east regions. 


these authors found evidence for a more extended emission 
with a much softer spectrum. 

On the basis of a new model for the synchrotron com¬ 
pact nebular emission, Sefako and De Jager (2003) esti¬ 
mated for the Vela plerion the acceptable ranges of the 
two main PWN parameters: the wind magnetization a 
(Kennel & Coroniti 1984), whose value is expected in the 
range 0.05-0.5, and the pair production multiplicity M, 
in the interval between 300 and 700. These authors used 
the power law spectrum derived from OSSE data and as¬ 
sumed that the Spectral Energy Distribution (SED) peaks 
in the MeV range. 

In this paper we present the results of the analysis of 
BeppoSAX data of the Vela Plerion in the energy range 
3-200 keV combined with the imaging and spectral capa¬ 
bilities of XMM-Newton. We derived information on the 
size of the source emission above 10 keV and obtained 
an accurate evaluation of the spectral parameters, useful 
to measure the total energy output from this plerion and 
for a more accurate comparison with models of nebular 
emission. 

2. BeppoSAX and XMM-Newton observations 
and data reduction 

The XMM-Newton observations were performed on 1 and 
2 December 2000 with an exposure of 38 and 57 ks, respec¬ 
tively. We reduced the EPIC data with the XMM Science 
Analysis System (SAS) version 6.0 and performed stan¬ 
dard screening of the EPIC data. The 2 December 2000 
observation suffered high background contamination and 
required a careful rejection of time intervals with soft pro¬ 
ton flares; to asses the degree of soft proton contamina¬ 
tion we used the diagnostic defined by De Luca & Molendi 


(2004). Because of this additional screening the effective 
exposure in the 2 December 2000 observation was reduced 
to 47 ks. In both observations MOSl was operated in 
Large Window mode, where the inner active region of the 
MOS array is the central square of 5'.5 side (see Fig. 1, 
left panel). MOS2 was working in Small Window mode 
with an active central square of lb 8 side only and was 
not included in the analysis of the wider central plerion 
emission. PN was working in Small Window mode but 
the target was not centred in the active square region and 
therefore PN data were not included in our analysis. 

Diffuse and asymmetric emission is visible in the im¬ 
age of Fig. 1 (left panel) beyond 8' mostly in South-East 
and South-West directions. Comparing the X-ray contour 
with the radio images recently obtained by Dodson et al. 
(2003) we note that the structure of the diffuse emission 
is correlated with that observed at radio frequencies. We 
extracted spectra from circular regions centred on the 
source with a radius of O'.5, I'.O, 1'.5 and 2'.0 (there¬ 
fore including the central source) and, to investigate bet¬ 
ter the spectral changes moving outward in the nebula, 
also from annular regions centred at the source with radii 
0'.5-T.0, T.0-T.5, l'.5-2'.0 and 8'-12' (outside the gap). 
We also extracted spectra from two circles of 2' radius 
in the South-East (RA: 08h35m52s, DEC: -45dl6m26s) 
and South-West (RA: 08h34m58s, DEC: -45dl7m04s) re¬ 
gions where, according to Fig. 1, diffuse emission is visible. 
Spectra and effective area corresponding to the two obser¬ 
vations were summed using the tasks mathpha and addarf 
from the HEASOFT 5.3.1 package. 

The spectrum of the background for the MOSl was 
extracted from a set of 10 public galactic plane observa¬ 
tions at Z=310°-312° from which we have removed the de¬ 
tected sources, and, for comparison, from closed datasets 
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(Marty et al. 2002) containing only the instrumental and 
particle components of the background. Results obtained 
using these two different background spectra were always 
compatible. Vignetting correction was properly applied to 
both source and background spectra. Spectral response 
matrices and effective area were calculated for each re¬ 
gion via the SAS tasks rmfgen and arfgen, respectively, 
according to the standard recipe for vignetting-corrected 
extended sources. 

The BeppoSAX. observation of the Vela pulsar region 
was performed on 18-20 November 1997 for a total net 
exposure of 84 ks in the MEGS and 38 ks in the PDS (for 
details on the BeppoSAX instrumentation see Boella et 
al. 1997). Standard procedures and selection criteria were 
applied to the data using the SAXDAS v.2.0.0 package. 
The image extracted from the MEGS (Fig. 1, right panel) 
presents regions of diffuse emission in the South-East and 
South-West directions just as in the XMM image. The 
extended structure of the source requires that a suitable 
auxiliary matrix must be created to correct for the re¬ 
distribution due to the instrument point spread function 
and for the vignetting relevant at distances larger than 
about 4' from the centre of the FOV (Fiore et al. 1999). 
A proper response matrix can be generated for the MEGS 
data using the SAXDAS effarea command (for details see 
Molendi 1998, D’Acri et al. 1998) and the average radial 
profile of the Vela plerion derived from the XMAi-Newton 
image interpolated with a spline function in the gap region 
and extrapolated with a constant up to 15'. We produced 
ad-hoc matrices for the circular collecting areas of radii 
equal to 4', 8 ', 12' and 15'. The poor angular resolution of 
the MEGS and the difficulty of producing a good response 
matrix for annular regions did not allow us to perform a 
spatial differential study of the X-ray emission, as done 
for MOS data. 

The spectrum of the background for the MEGS was 
extracted from archive blank fields. 

PDS data were collected using the standard rock¬ 
ing observation mode with two detectors pointing at the 
source and the other two in offset directions to measure 
the background simultaneously. The net count rate in the 
whole energy band was 1.7 c/s. The PDS field of view, 
1 . 3 ° FWHM, is large enough to include the entire ple¬ 
rion emission. It is, however, much narrower than that of 
OSSE ( 3°.8 X 11°.4 FWHM). We also checked that PDS 
field of view does not contain bright sources that could 
contaminate the Vela spectrum. 

Events with an energy higher than 3 keV were selected 
for the analysis to exclude the soft thermal contribution 
from the pulsar (Pavlov et al. 2001b) and from the Vela 
shell (Bocchino et al. 1999). Effective ranges were 3-10 
keV for MEGS and MOS and 15-200 keV for PDS. For 
the same reason we did not consider LEGS data in the 
analysis. Spectral fits were first evaluated separately for 
the MOS, MEGS and PDS data and after, combined fits 
of MOS-PDS and MEGS-PDS data sets were considered. 

A relevant point we stress when working with spectra 
obtained with several instruments is that a proper evalu¬ 


ation of the inter-calibration factors is required. For the 
BeppoSAX NFIs, the accurate ground and in-flight cali¬ 
brations were used to establish that the admissible ranges 
for the factor between MEGS and PDS for point sources 
(/mp) is 0.77 < /mp <0.93, reduced to 0.86±0.03 for 
sources with a PDS count rate higher than 2 ct/s (Fiore 
et al. 1999). For extended sources, as for the Vela plerion, 
we expect that the best fit estimate of /mp would result 
in agreement with the above value only when the emission 
observed by the PDS is fully selected in the MEGS image. 
Moreover, MOS-PDS intercalibration factor is expected 
in a similar range because the intercalibration analysis be¬ 
tween MOS and MEGS shows a good agreement between 
the two instruments (Kirsh 2004). Values of the intercal¬ 
ibration factors well above the expected range can be an 
indication that MEGS flux, for a given collecting area, in¬ 
cludes only a fraction of the signal detected by the PDS. 
Thus it is possible to use this parameter to obtain infor¬ 
mation on the angular size of the source responsible for 
the hard X-ray (>15 keV) emission. 

The uncertainties reported in the following are at 1 
standard deviation for one interesting parameter. 

3. Spectral analysis 

3.1. Single instrument spectral analysis 

Spectra of each instrument were fitted with a single power- 
law model. A plot of the MOS and MEGS photon index 
as a function of the extraction radius for all the circular 
regions centred at the source is shown in Fig. 2. The value 
of r changes with the extraction radius: for a radius <2' it 
is around 1.55, while for higher radii it increases to 1.67, in 
agreement, within the errors, with the value reported by 
Willmore et al. (1992), indicating a quite mild softening 
of the spectrum in the outer region of the PWN. Figure 
3 shows the fluxes in the 3-10 keV band as a function 
of the spectra extraction radius. The flux increase within 
4' quantifies the amount of emission originating from the 
inner region. 

The spatially differential spectral analysis confirms 
that the inner region of the plerion has a flatter spec¬ 
trum than the outer one. The average spectral index of 
the three annuli is 1.60±0.02 (the third annulus differs 
from this value only by two standard deviations) to be 
compared with 1.90 found at distances greater than 8 '. 
The small difference (at about two standard deviations) 
between the mean index of the three inner annuli and that 
found for the circular region of O'.5 radius is likely due to 
the contribution of the inner part of the PWN which is 
about 40% of the total as derived from the 3-10 keV flux 
ratio. 

The spectral index variation in the inner three MOS 
annuli is in agreement with the Ghandra spectral index 
map presented by Kargaltsev & Pavlov (2003). 

Fitting the PDS spectrum in the energy range 15—200 
keV gives a photon index equal to 2.00±0.05 with the Xr 
of 0.73 for 16 d.o.f. and a flux of (15.7±0.4)xl0“^^ erg 
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Fig. 2. Panel (a): the power law photon index in the 
3—10 keV band of MOS (stars) and MEGS (filled circles) 
spectra at varying extraction radius. Panel (b): the power 
law photon indices in the 3—10 keV band of MOS spectra 
extracted from the annuli listed in Table 1. 
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Fig. 3. Vela PWN flux in the 3-10 keV band at varying 


extraction radius. 


cm“^ s“^. This photon index, significantly higher than 
those found in the lower energy range, indicates that a 
spectral steepening occurs at higher energies. Moreover, it 
is higher than the value obtained with OSSE (De Jager et 
al. 1996), although statistically compatible because of the 
large uncertainty in their estimate. A summary of the best 
fit results for the various extraction radii of the imaging 
instruments and PDS are given in Table 1. 

We also investigated the spectral properties in the 
outer part of the plerion. The XMM spectra extracted 
from two circular regions having a diameter of 4' in the 
South-East and South-West were also fitted with a simple 
power-law model. We obtained P = 1.9 ± 0.1, K = (5.0 ± 
1.7) X 10“^ photons cm“^ s“^ keV“^ (Xr = 0.94, 181 d.o.f) 
in the SE region and P = 1.9 ±0.1, K = (5.5 ± 1.5) x 10“"^ 
photons cm“^ s“^ keV“^ (Xr = 1-01, 186 d.o.f) in the SW 
region, respectively. The 3-10 keV band fiuxes per unit 


Fig. 4. Best fit values of /mp in the combined log- 
parabolic fit with the PDS at varying extraction radius. 
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Fig. 5. Best fit of MEGS and PDS data with a broken 
power law model (15' extraction radius for the MEGS 
data). See Table 2 for the best fit parameter values. 


3.2. Combined instrument spectral analysis 

The finding that the 15-200 keV X-ray photon index is 
significantly larger than those measured below 10 keV in¬ 
dicates that a single power law is not able to fit the broad 
band spectrum of the Vela PWN. We therefore performed 
a spectral analysis combining the PDS data with those of 
MOS/MEGS. The single power-law model was rejected be¬ 
cause the Xr values of the combined fits were in the range 
1.34(341 d.o.f.) - 1.84(83 d.o.f.). Moreover, PDS residuals 
show very large deviations with Xr always higher than 4 
(15 d.o.f). 
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Table 1. Best fit parameters for a single power law of MOS and MEGS. The values given in the first section are for 
the spatially integrated analysis, while those of the second section are for the MOS annular regions. 



Radius 

p(a) 

^(a,b) 

p{a,c) 

Xr (d.o.f.) 

MOS 

0'.5 

1.50±0.02 

0.39±0.01 

1.785±0.009 

1.16 (253) 

MOS 

T 

1.55±0.01 

0.95±0.02 

3.98±0.01 

1.08 (324) 

MOS 

T.5 

1.56±0.01 

1.09±0.02 

4.49±0.01 

1.11 (331) 

MOS 

2' 

1.55±0.01 

1.18±0.02 

4.95±0.01 

1.07 (337) 

MEGS 

4' 

1.61±0.02 

1.50±0.04 

5.66 ±0.03 

0.96 ( 66) 

MEGS 

8' 

1.67±0.02 

1.88±0.05 

6.39T0.03 

0.93 ( 66) 

MEGS 

12' 

1.66±0.02 

2.17±0.05 

7.50T0.03 

1.02 ( 66) 

MEGS 

15' 

1.66±0.01 

2.62±0.09 

9.06T0.04 

1.02 ( 66) 

PDS 

FOV 

2.00±0.05 

6.19±0.95(‘^> 

15.9±0.4(®> 

0.73 ( 16) 

MOS 

0'.5-T 

1.60T0.02 

0.56T0.02 

2.150T0.009 

1.05 (209) 

MOS 

T-l'.5 

1.61±0.03 

0.14T0.01 

0.528T0.005 

0.88 ( 72) 

MOS 

T.5-2' 

1.50±0.05 

0.08±0.01 

0.361T0.004 

0.77 ( 49) 

MOS 

8'-12' 

1.90±0.06 

0.51±0.05 

1.165±0.02 

1.13 (231) 


Errors correspond to 1 standard deviation for one interesting parameter. 

Power law normalisation K is expressed in units of 10”^ photons cm“^ s~^ keV~^ at 1 keV. 

Flux is in units of 10”^'^erg cm“^ s“^. It is relative to the 3-lOkeV band for the MOS/MECS and to the 20-100keV band 
for the PDS. The reported errors are statistical. Systematical errors amount to 5% for MOS flux values and 10% for the 
MEGS. 

The PDS normalisation is obtained dividing the single power law fit result (5.2±0.8) by the intercalibration factor /mp = 0.84 
obtained in the combined fit with the MEGS spectrum extracted from the circle of 15' radius. 

The PDS flux value is that obtained from the single power law fit (13.39±0.34) divided by the intercalibration factor 
/mp = 0.84. 


Two other models were then considered: a broken 
power law: 

F{E) = K {E/1 keV)"^i for E < Eb 

F(E) = (E/lkeV)-^" ior E > Eb (1) 

and a continuous steepening law, characterised by a linear 
dependence of the spectral slope upon the logarithm of 
energy: 

F{E) = K{E/1 keV)“(“+^ (B/ikeV)) (2) 

The fits of the broken power law (Eq. 1) to the spec¬ 
tra for different extraction radii gave acceptable Xr for 
all regions but showed a strong correlation between the 
intercalibration factor /mp and the break energy Eb- To 
extract reliable values from this model it is necessary to 
fix one of these two parameters. 

The log-parabolic law (Eq. 2) gave acceptable fits for 
all extraction radii without evidence of a correlation be¬ 
tween the parameters. The best fit values of /mp are out¬ 
side the expected range with extraction radii lower than 
15' as shown in Fig. 4. This result indicates that to match 
the PDS signal to the 3-10 keV flux we must extract events 
in a MEGS region with a radius of about 15'. 

To further investigate the extension of the region re¬ 
sponsible for the emission above 15 keV we fixed /mp 
to 0.93 (the highest acceptable value) and, fitting the 12' 
spectrum with Eq. (2), we obtained an acceptable Xr of 
0.98 (83 d.o.f.) and an acceptable Xr relative to PDS data 
only of 1.09 (15 d.o.f). However, the same procedure on the 


8 ' spectrum gave a marginally acceptable fit for the two 
datasets {Xr = 1-3, 83 d.o.f.) but not acceptable for PDS 
data only {x/ = 2.90, 15 d.o.f.). Thus, we are strongly 
confident that the extension of the hard X-ray emitting 
region is at least 12' from the pulsar. 


Table 2. Best-fit model parameters of the combined fit 
MEGS-PDS relative to the MEGS extraction radius of 15' 


Log-parabolic Model 
a 1.36±0.05 

(5 0.21±0.03 

K (x 10"^) 2.08±0.10 

/mp 0.84±0.03 

xl (d.o.f.) 0.76 (82) 


Broken power-law Model 


El 

Ea 

Eb 

K (x 10"^) 
Imp 

Xr (d.o.f.) 


1.66±0.01 
2.01±0.05 
12.5±1.5 
2.62±0.06 
0.84 (fixed) 
0.96 (82) 


We also fitted the MEGS spectrum extracted from a 
15' radius and the PDS spectrum with the broken power- 
law model, fixing the intercalibration factor to the value 
obtained in the log-parabolic fit. The fit gave an accept¬ 
able Xr- The spectrum and the best fit model are shown 
in Fig. 5, and the parameters’ values are given in Table 2. 
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Note that the hard branch of the spectrum is fully consis¬ 
tent with the fit of PDS as single instrument (see Table 

!)■ 

We then verified if the models used to describe the 
spectral distribution of the Vela PWN in the 3-200 keV 
were consistent with the EGRET upper limits reported 
by Kanbach et al. (1994) for the unpulsed component ob¬ 
served from the region of PSR 0833-45. We extrapolated 
the two best fit spectra of the combined MECS-PDS data 
sets (15' extraction radius for the MEGS) in the EGRET 
range and the results are shown in Fig. 6 . While the log- 
parabolic spectrum is much lower than the 7 -ray upper 
limits, that derived from the broken power law is only 
marginally compatible with them. It is possible, however, 
that the plerion spectrum will show another break and/or 
a cut-off at energies in the range between 200 keV and 
~10 MeV. In the same figure we reported also the upper 
limits of Mignani et al. (1989) for the optical nebular flux 
which are in satisfactory agreement with the low frequency 
extrapolation of the inner plerion spectrum. 

4. Discussion 

The broadband X-ray spectrum of the Vela PWN de¬ 
rived from the combined analysis of BeppoSAX and XMM- 
Newton data, together with the information from high res¬ 
olution images of Ghandra (Pavlov et al. 2001a,b, 2003), 
provided a more complete detailed picture of the spec¬ 
trum and the structure of this source useful for the un¬ 
derstanding of the main physical processes responsible for 
the energy transfer and radiation. 

We can distinguish two main regions where X rays are 
emitted: a central source of high brightness with an angu¬ 
lar diameter of 3'-4' surrounded by a much fainter region 
extending mainly in the SE and SW sector up to distances 
of about 12'-15', corresponding to a distance from the 
pulsar of about 1.2 pc. This picture has a good correspon¬ 
dence to the radio images recently obtained by Dodson 
et al. (2003). These authors report the presence of an ex¬ 
tended and highly polarised non-thermal emission out to 
their held limit of about 6 ' from the pulsar. The spatial 
structure is mainly elongated in the NE-SW direction, but 
towards the NE it is limited by a quite sharp boundary, not 
detected in the other direction. Furthermore, the images 
of Dodson et al. (2003) suggest the presence of another ex¬ 
tended component in the SE direction at distances larger 
than ~ 6 '. Likely, these two outer components are the ra¬ 
dio counterparts of the extended X-ray emission detected 
in the MOS and MEGS images. 

Our spectral analysis conhrms that the central region 
has a harder spectrum than the outer one, with a photon 
index of about 1.6, but possibly close to 1.5 as indicated by 
the MOS spectra with the smallest extraction radius. Such 
a hard spectrum was also observed in the high resolution 
X-ray images obtained by Pavlov et al. (2003) with the 
Ghandra observatory. The photon index from MOS data 
in the outer region is indeed close to 2 . 


The spectrum at energies higher than ~15 keV is well 
fitted by a single power law with a photon index 2. A 
possible interpretation of these results is that the emission 
from the extended region is described by such a spectrum, 
while that of the central PWN must show a break at ~12 
keV with a photon index change of about 0.5. 

This behaviour does not confirm the claim of De Jager 
et al. (1996) about the spectrum of the central compact 
nebula extending unbroken up to 400 keV. 

In Fig. 6 we compare our broad band results with the 
hard X-ray SEDs observed by Birmingham 2 (Willmore 
et al. 1992) and OSSE-GGRO experiments (De Jager et 
al. 1996). Our spectrum is well above both these data, 
however, the discrepancy with the data of Willmore et al. 
(1992) can be explained by the different source extension 
taken into account. The data of the Birmingham 2 ex¬ 
periment were obtained modelling the coded mask image 
with a central compact source plus an extended emission. 
The SED of Fig. 6 refers to the central source only and 
when scaling our SED to the flux measured by the MEGS 
with an extraction radius of 4' (see Table 1) we found 
a good agreement. At variance, the discrepancy with the 
spectrum derived from the OSSE-GGRO observation (De 
Jager et al. 1996) is not simply explained by the source ex¬ 
tension and can be due to some systematic effect possibly 
due to the much wider field of view of this experiment. In 
particular, when these data are extrapolated from 60 keV 
to energies of ~1 keV (Sefako and De Jager 2003), the 
source X-ray luminosity can be largely underestimated. 
We computed the isotropic X-ray luminosity of the Vela 
plerion by integrating between 0.1 and 200 keV the bro¬ 
ken power law spectrum taken over the largest radius and 
obtained Lx = 5.5 x 10^^ erg s“^. This luminosity is more 
than a factor of 2 higher than that derived from the best 
fit OSSE spectrum in the same range. Note that about 
45% of this luminosity is emitted at energies lower than 
Eb. 

A comparison of the PDS flux of the Grab in the 20- 
200 keV range and the corresponding OSSE flux based on 
Ulmer et al. (1995) also shows a discrepancy of the order 
of 30%. 

Pavlov et al. (2001a) estimated that the 1-10 keV lu¬ 
minosity of the PWN within about 1' from the pulsar is 
6 x 10 ^^ erg s“^, in a good agreement with that derived 
by us from the MOS data for a similar extraction radius 
(Fig. 3). Such a luminosity corresponds to about 1/4 of 
the emission from the whole source and implies that the 
outer nebula emits a relevant fraction of the total X-ray 
power. This more precise value can be useful to improve 
the estimates of the magnetization parameter a and the 
pair multiplicity M of the Vela plerion. We did not develop 
a complex physical model of the source for this purpose, 
which is beyond the scope of the present paper. However, 
one can reasonably expect that the greater luminosity ob¬ 
tained by us would imply a higher number of particles and 
consequently a higher value of M. If confirmed by numer¬ 
ical models, like that by Sefako and De Jager (2003), this 
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Fig. 6. The optical to X-ray spectral energy distribution of the Vela PWN compared with other literature results. 
Our SEDs are the broken power law (solid and short-dashed lines) and the log-parabola (dotted line) to which 
are superposed the deconvolved spectral data. The thick solid line is the spectrum measured by the Birmingham 
2 experiment (Willmore et al. 1992) and the dotted lines represent the broken power law SED rescaled to the flux 
measured by the MEGS with an extraction radius of 4'. The thick long dashed line is the broad band spectrum of 
the compact nebula based on OSSE data (De Jager 1996) and used by by Sefako and De Jager (2003). Optical upper 
limits (empty triangles) are taken from Mignani et al. (2003) extrapolated to an area of 10"^ arcsec^ to cover all the 
inner part of the plerion. The triangles in the 7 -ray band represent EGRET upper limits from Kanbach et al. (1994). 


finding can be the reason for the too low value found by 
these authors. 

Furthermore, we found that the X-ray emission from 
the plerion Lx is a factor of the order of ^ 10 “^ of the 
pulsar spin-down luminosity E = 6.7 x 10^® erg s“^, one 
order of magnitude higher than the value given by Pavlov 
et al. ( 2001 a). 

It is possible to show, using a simple dimensional ap¬ 
proach, that the nebular extension measured in the MEGS 
images and used in our analysis (i.e. a radius 12'-15', 
approximately corresponding to 1 . 0 - 1 .3 pc) is compati¬ 
ble with the radiative lifetime of high energy electrons 
tsy emitting hard X rays. According to a general plerion 
model we assume that electrons are injected at the pulsar 
wind shock, which for Vela is at an angular distance from 
the pulsar of 53" (Helfand et al. 2001), and then diffuse 
throughout the plerion. From the synchrotron theory the 
radiative lifetime of the electrons which emit photons of 
energy is: 


where is the average transverse component of the mag¬ 
netic field. The typical time scale to escape the nebula of 
radius R (in pc units) can be written as 


te — R/vd = 1.03 X 10® a R 


( 4 ) 


where 1 /a is the outward diffusion velocity in units of the 
speed of light {vd = c/a); for example, if Vd is equal to the 
sound velocity in a relativistic plasma then a = (De 
Jager et al. 1996). The two above equations can be used 
to estimate the highest energy of the photons emitted by 
electrons able to reach the outer boundary of the plerion 
(i.e. te — tsy) and then we obtain: 


, _ 435 / R \ Y Bx 

“ a2 Vlpcy Vl0”^G 


-3 


keV 


( 5 ) 


The typical nebular magnetic field is estimated to be 
about 10““^ G (Sefako & De Jager 2003, Pavlov et al. 2003) 
and therefore, for any reasonable choice of a, photons of 
energy of ^ 100 keV can be efficiently radiated up to dis¬ 
tances of about 1 pc from the wind shock region with¬ 
out need for a further acceleration. Note that these rela¬ 
tions imply a maximum Lorentz factor for the electrons of 
~3xl0®, as also estimated in the above quoted papers. 

An interesting consequence of Eq. (5) is that photons 
with energy of ^ 10 MeV should be emitted from regions 
much closer to the wind shock and therefore the outer 
nebula spectrum should show a further steepening at en¬ 
ergies higher than a few hundred keV. The occurrence of 
such a high energy break could be detected by the study 
of the DG component of the Vela pulsar with the present 


s 
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(INTEGRAL) and the next generation of high sensitivity 
low-energy 7 ray telescopes. 
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